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I. INTRODUCTION 


A. EVOLUTION OF EARTHQUAKE PREDICTION 


The problem involving searching for a reliable precursor 
to earthquakes is being given much attention throughout the 
world. With the reality that within minutes an earthquake 
can trigger an incredibly disasterous act of nature leaving 
death tolls in the tens of thousands and property damage in 
the biliions anywhere in the world. The desire for 
prediction capabilities has evolved into an intense investi- 
gative scientific undertaking. Until recently, seismic 
activity data was collected cnly to establish where an 
earthguake has the greatest potential to occur. Hundreds of 
potentially destructive seismic shocks escape unnoticed 
because they occur in remote, sparsely populated areas of 
the earth. 

The present state of the art in earthquake prediction 
may be compared to that of weather forecasting some 50 years 
ago. Seismological data acquired must be continuous and 
detailed from every area of the world that is subject to 
seismic disturbances. Various experimental methods for 
prediction of seismic activity are widely used such as elec- 
tric and magnetic field, tilt, creep, and strain motion 
measurements to establish a practical data base [Ref. 1]. 

Evidence of radio (RF) emission associated with precur- 
sors to earthquake activity is very sparse. Electromagnetic 
radiation associated with earthquakes appears impulsively 
[Ref. 2], which suggests that these electromagnetic waves 
are radiated in a wide range of radio frequencies. The char- 
acteristics of the electromagnetic radiation depend greatly, 


however, on such factors as the magnitude of the earthquake, 
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the depth of focus, the geology of the earthguake region, 
the distance between epicenter and point of observation, the 
condition of Lonospheric propagation, and others. This is 
possibly the reason why the intensity of the electromagnetic 
radiation is different for different earthguakes, with 
values varying from tens or hundreds of microvolts per meter 
up to tens of hundreds of volts per meter. 

Measurements are presently being made to explore the 
possibility of the existence of electro-magnetic fields. 
However, there presently appears to be no directed research, 
or a plausible explanation for, reported radio frequency 
emission interference coinciding with an increase in earth- 
guake activity. Derr noted an okservation of noise on the FM 
broadcast band (88-108 MHz) cccuring Simultaneously with 
earthguake lights [Ref. 3]. The reception of earth-based 18 
MHZ radio noise by JW. Warwick of the High Altitude 
Observatory, Boulder, Cow. at the time of the famous Chilean 
earthquake still has been unanswered [Ref. 4]. Research 
personnel involved in RF interference study, in intervie 
citizen band and amateur radio operators in the Hollister, 
Ca. area, noted a consistent report of "increase in back- 
ground noise level preceding earthguake activity" in the 
freguency range of CB (27 MHz) and six and two meter amateur 
bands (50-54 MHz and 144-148 MHz). Reports typically cited a 
noticeable (10 = dB) increase in background noise which 
preceded an earthguake by 12 to 24 hours, returning to 
normal just before the tremors were felt (one hour to 15 
Minutes). Observations of the electrostatic field (field 
mills) have reported that diurnal variations at selected 
Sites are accompanied by "high frequency, spiky noise" that 
apparently have no explanation. Various theories have been 
approached on what is actually occuring at the fissure 
boundary, whether impact of dissimilar materials oor shock 


induced polarization have developed the necessary charge 
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density to produce the electric field and subsequently the 
specific surface area behaves as a constant current 


generator. 


Be. FOCUS OF STUDY 


The design, construction, and testing of a freestanding 
radio frequency (RF) noise measurement and recording system 
in the 30 MHz and 150 MHz range is presented. Placement of 
the system along the San Andreas fault to permit the estab- 
lishment of a correlation between increased background RF 
noise to earthquake fault activity for potential as an 


earthquake prediction tool is described. 


1s 


II. BACKGROUND 


A. EARTHQUAKE GROUND MOTION 


Equally important to predicting the occurence of an 
earthguake is determining which of the many ways the ground 
is likely to shake during the earthquake, how strong the 
Shaking will be, and how long it will last. Knowledge of 
ground motion can make it possikle to understand the charac- 
teristics of the earthquake source. Most measurements made 
during actual earthquakes at seismological stations are some 
distance from the source. The localized source of an earth- 
quake was not observed for its significance until the San 
Francisco earthquake of 1906. 

It was recognized that slippage along a fault in the 
earth's crust caused the earthguake. The classic study of 
Ha Geyer. Reid of John's Hopkins University revealed that 
tocks parallel to the fault had been strained or sneared. 
From the basis of his geodetic observations eminated the 
elastic rebound theory of earthquakes which asserts’ that 
rocks are elastic, and mechanical energy can be stored in 
them just as it is stored in a compressed spring {Ref. 5]. 

When two Llocks forming the opposite sides of a fault 
move by a small amount, the motion elastically strains the 
rocks near the fault. The frictional bond fails at its 
weakest point when the stress becomes large enough. That 
point of initial rupture, called the hypocenter, may be near 
the surface or deep below it. The seismic waves radiate from 
the hypocenter in all directions producing the earthquakes. 
The point on the surface of the earth above the hypocenter 
is the epicenter of the earthquake. The rocks rebound over 


an interval of minutes, days or even years. The seismic 
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energy radiated over ground at any one time may be quite 
small. This slow process is referred to as a seismic slip or 
creep. 

Most earthguakes are generated in zones where the huge 
plates of the lithosphere, which make up the outer layer of 
the earth's surface, as in Figure 2.1 and Figure 2.2 
[Ref. 6], are shearing past each other. A conceptual model 
of the faulting process, as in Figure 2.3 and Figure 2.4, 
were largely developed from investigations of earthquakes 
along the San Andreas fault. Over the past 137 years, more 
than 20 severe earthguakes have been recorded in the state 
of California, mainly along the San Andreas fault. The San 
Andreas fault which extends frem the Gulf of California in 
the South to Cape Mendocino varies in width from about a 
meter to 1.6 KM, and plunges almost vertically into the 
earth for 15 to 20 KM. These are some of the reasons why the 
San Andreas qualifies for conducting an intensive, longtern 
RF noise measurement investigation to show anomalous 
increases in the intensity of electromagnetic radiaticn 


precursor to earthquakes. 


Be SOURCES OF ELECTROMAGNETIC RADIATION IN EARTHQUAKES 


The emissions of RF are theorized to be caused by the 
increasing of mechanical pressure at the fault causing shear 
stresses that convert the insulating rock to a conductor {or 
semiconductor) by reducing the energy gap between the 
valence and conduction bands. During the stress increase, 
microfractures cause shock waves that accelerate the elec- 
trons in the conducting bands to radiate both acoustic and 
electromagnetic energy. If the mechanical forces increase to 
the yield point, the rock fails in brittle fracture and 
causes a sudden decrease in the shear forces. The mechanical 


motion is described by a stick-slip motion that has stable 
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Figure 2.2 Obligue View Across the San Andreas Fault Zone 
Typical Right-Lateral Offset of Strean 
Offset is About 450 Feet 
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Figure 2.4 Double Couple Source Model of Fault Slippage 
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1. Insulators and Conductors 


The radio emission theory restS onthe principle 
that an insulator can be mechanically stressed to become a 
conductor or semiconductor with subsequent charge migration 
from the conducting rocky material to the insulating rock 
under the equivalent electric field generated by a_ shock 
wave during brittle fracture of the rock. Experimental 
evidence from the Ballistic Research Laboratory, Aberdeen 
Proving Ground, Maryland showed that shock-induced polariza- 
tion is an increasing function of pressure [Ref. 7]. 

At low pressures shock-induced polarization (Ps) 
becomes small but 1s easily measured down to approximately 
25 Kbar. Shock-induced polarization of both plexiglas and 
polystyrene was observed at lower pressures, but the noise 
level associated with the experiments prevented accurate 
measurement. At high pressures data indicated transitions Ps 
began a more rapid increase, suggesting transitions toa 
State where the plexiglas and polystyrene were more easily 
polarized by the shock wave. 

The dielectric constant of the shocked mediun, K'*, 
was also an increaSing function of pressure. The values for 
K' were characterized by appreciable scatter, but indicated 
a gradual increase up to the transition pressure where a 
more rapid increase was noted. The relaxation time, “) , was 
also characterized by appreciable scatter, but remained 
between 0.5 and 1.5 usec up to the transition pressure, 
where a rapid decrease began. If it is assumed that polari- 
zation is caused by the motion of atoms or atomic groups, 
then the sudden decrease in “T suggest that motion cccurs 
more freely in the high pressure phase, and that the 


transitions are to a more fluid state. Greater freedom of 
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motion would result in a larger polarization of the material 
by both mechanical and electrical forces in aggreement with 


the observed increase in both Ps and K*t at the phase 
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Figure 2.5 Shock-Induced Polarization Signals 
of Plexiglas Below Transition Pressure 


transition. Figures 2.5, 2.6, 2./, and 2.8 show that over 
the entire pressure range of the measurements, the polariza- 
tion signals can be described by a physical model involving 
three variables, Ps, K'*, and ‘Pf . It should be noted that 


the computed curves are in better agreement with measured 
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Pigure 2.6 Shock-Induced Polarization Signals 
of Plexiglas Above Transition Pressure 


Signais at low pressures than at high pressures where the 
ends of the measured signals are more distorted by curvature 


and obliquity of the shock front. 
2. Contact Blectripacaien 


Electrification of insulators and conductors ~b9 
impact are described in experiments conducted by Schnurmann 
and Warlow-Davies of two metals separated by an insulating 


orie: Slow sliding can proceed by a series of jerks 
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Figure 2.7 Shock-Induced Polarization Signals 
of Polystyrene Above and Below Transition Pressure 


exhibiting two characteristics, that of relaxation oscilla- 
00 . The first is a fixed amplitude and the second is a 
period determined by a relaxation time for the building up 
of a limiting value. The experiments described the amplitude 
of the jerks was observed to decrease when the dielectric 
breakdown strength of the material used aS a lubricant was 
diminished, or when, with the same lubricant, the velocity 
of Sliding was increased. Observations of contact 


electrification show that an electrostatic force resisting 


2S 








© - PREXIGLAS 
O - POLYSTYRENE 


X\ 


| 
! 
in : 

oe ) 
/ 





ee ee ae 


a 


US 156 U es 300 
PRESSURE = Kbor 





ee | 


Figure 2.8 Shock-Indaced Polarization wake of Polystyrene 
and Plexiglas as a Punction of Shock Pressure 


sliding must exist if the boundary layer consists of dielec- 
tric or semiconducting matter and this electrostatic compo- 
nent is an appreciable part of the force of sliding friction 
[Ref. 8]. 

Active fault lines move by relaxation oscillations 
and can be described by the nonlinear stick-slip notion that 
is generated when the static force of friction is larger 
than the kinetic frictional value. During the increasing 
stress, the motion of the fault proceeds uniformly by creep 
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and the large stress converts the insulating rock to a 
conductor or semiconductor. On reaching the brittle fracture 
stress, the material unloads and the stress is reduced and 
results in a retreating shock wave inthe material. The 
shock wave has a conductor ahead of it with free electrons 
that can fall into the insulatcr conduction band behind the 
retreating shock front. The charges can accumulate by succe- 
Sive creep-slip moticns to cause the electronic charge to 
produce an electrical field beyond the breakdown value in 
the air or rock surface, cauSing an arc with the attendant 


Tadio emissions. 


3. Fracture of Quartz-bearing Rocks 


ae RE Sa a Se a 


The geological research group at the Los Alamos 
Scientific Laboratory were believed to be the first to 
detect radio frequencies when gquartz-bearing rocks and other 
hard piezoelectric materials were fractured in tJlaboratory 
experiments. Using a standard RF coil with a ferrite core 
placed about 5 to 50 CM from the sample and a transient 
recorder it was found that small cracks in quartz and tour- 
mMaline single crystals produce exponentially decaying RF 
transients with a time constant of about 1075 seconds as 
shown in Figure 2.9 . Rock samples were fractured by 
applying uniaxial load between two parallel steel anvils or 
by pressing a small steel ball against a flat surface of the 
Sample. The signal developed on the coil was amplified by a 
wideband amplifier and fed into. one channel of the transient 
recorder. All measurements were performed under ambient 
conditions. Also the samples were not treated in any way 
except for cutting or coring [Ref. 9]. 

The fact that the fracture related EM emissicn was 
observed only in quartz and tourmaline crystals and in 
quartz-bearing rocks rules out triboelectric effects 
(frictional electricity) as the source of this emission and 


indicates that it is related to piezoelectricity. The most 
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Fig. 2. Spectral content of fracture-related 
EM emission from (mean grain size shown in 
brackets): a. coarse-yrain granite (15 mn); b. 
granodiorite (2 mn); c. fine-grain granite 
(from Westerly, R. I., 0.5 mm); d. fine-grain 
quartzite (0.1 me)ijes de7zax ete ec emem 
tourmaline single crystal cracked by arplying 
pressure along the c axis. The curves show 
averages of 10 spectra for each sample. 
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Figubewz Spectral Content of Fracture-Related EM Emission 
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probable mechanism for this emission 1s the rapid drop in 
the piezoelectric field accompanying the sudden release of 
stress when fracture occurs aS Shown in Figure 2.9 . This 
hypothesis is supported by the observation that the ampli- 
tude of the resulting RF signal is independent of MJload, 
since regardless of the stress distribution within the 
sample fracture occurs only when the local stress exceeds a 
specific value (much higher than the average stress in the 
experiments) and therefore drop in stress is essentially 
independent of load. 

Figure 2.10 compares the spectral content of frac- 
ture related EM emission from several samples, obtained by 
averaging 10 different spectra for each sample. It is clear 
that a shift to higher frequency occurs as the grain size of 
guartz-bearing rocks decrease. 

Detection and analysis of fracture-related RF 
Signals was very difficult at distances greater than 50 cn. 
With regard to RF emissions outside of the laboratory, no 
conclusions could be made from the study concerning whether 
a substantial increase in the size of the fracture results 
also in a corresponding increase in the emitted RF energy, 


thereby permitting detection at larger distances. 


4. Porous Rock 


Electrokinetic phenomera induced by ground water 
flow associated with earthguakes are believed to provide a 
possible means of earthquake prediction. Detectable varia- 
tions of earth currents, electric potential and geomagnetic 
field may be caused by diffusicn of fluid into a dilatant 
(wider or larger) focal region. The electrokinetic phenomena 
controlled by the processes of mineralized water migration 
through microcapillary rock channels are an independent 


source of charge separation and relaxation. The charge 


24 


Figure 2.10 


Fig. 1. Examples of EM transients accompanying 
fracture of piezoelectric specimens (the digi- 
tal sampling interval is given in brackets): a. 
small crack in a quartz singie crystal (50 ns); 
b. small crack in a tourmaline single crystal 
(10 ns); c. failure of a specimen of fame prain 
sandstone (200 ns); d. failure of a quartzite 
specimca, showing also the accompanying crop in 
load (2000 ns, see text). Note the different 
time scale for each transient. 
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separation is due to the selective absorbtion of ions fron 
the solutions on the fracture surfaces and the formation of 


double electric layers. 


J rid-Liquid Interface 


' 
x. = 


MO 





dee 


Yy 





Velocity Protile 


Ee potential 
Electric 


> 
Potential 


“ 
yy 
wae 








WW 
¥, 


= 


Vf 


Yt1. 


> Capillary 


+> 


ie 


. 


Electric Double Layer 


Neen ee ee aeons ieee nema aueianseecaenneemms ed 


Figure 2.11 Electric Double Layer and 
Velocity Profile in a Capillary 


The double layer is made up of a layer of ions which are 
firmly held to the solid, anda more diffuse mobile layer 
extending into the liguid phase. The resultant charge of the 
diffuse layer is equal in magnitude but of opposite sign to 
that of the firmly held layer. A simplified schematic 
diagram of the structure of the double layer and profiles of 
the electric potential and velocity of liquid in a capillary 
is shown in Figure 2.11 . Because of the electrical charges, 


there is a difference of electrical potential between the 


J, 


solid-liguid interface and the L-tulk of the liguid. This is 
called the electrokinetic potential or zeta potential. Waen 
an electromotive force (EMF) is applied to an electric 
double layer, there will be a displacement of the oppositely 
charged layers relative to one another. If the solid cannot 
move, the fixed part of the dcuble layer also is unable to 
move; the application of an EMF will result in movement of 
the diffuse layered ions and liguid [Ref. 10]. 


Applied EMF and viscous force is given by: 


= (eS /{)) grad E (2.19 


where j 1s the fluid flow flux (volume) through a capillary 
per unit time and per unit area, € isi the zeta potential, 
is the viscosity of the fluid, and € is the dieclectmee 
constant of the fluid. 

The streaming potential (E) is the production of a 
potential difference when a liguid is forced through a 


porous medium or a capillary tuke. It is given by 
Jraq.) =i (EX (Sg grad FE (2.29 


where and P are the electrical conductivity and the pres- 
sure of the fluid, respectively. 

Since all electrokinetic processes are irreversible, 
the thermodynamics of an irreversible process can be 
applied. The most general relations between the electric 
current 1 and the fluid flow j, and the forces of grad E and 


grad P are 
-l..= Lil’ gtad” b+ stiZeegeaae (2 a 


-j = £21 grad Ey ote i227 ora (2.4) 
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where Lij are constants. Explicit expressions of Lij's fora 


capillary tube are: 


L11= @ (2.5) 


L12 


bgi=-€5 a (2.6) 
L22 = k* // (en 


where k* is the specific permeakility of the capillary tube. 


For porous media k is given by 
= gk* (2.8) 


where d is the porosity of the porous media. The basic equa- 


tions of fluid flow in the porous medium are 


-I > & grad E- (PE\/N I grad P (29) 


“3 =- (P6979) grad B + (k /)) grad P (2. 10) 
It is assumed that the electrical conductivity of a matrix 
substance is much smaller than $0. 

The maximum magnitude of the current density associ- 


ated with the water flow is obtained by setting grad E = 0 : 
-I = (PES 79 ) grad P (2.11) 


This shows the electric current iS proportional to the pore 
pressure gradient (the pore pressure is assumed to change by 
10 - 1000 bars in the dilatant focal region preceding or 
following the earthguake). If the mean gradient of pore 


pressure is approximated to be 


grad P = 1= 100 bar/Kn (2.12) 
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the electric Current S@nauveowas 
iI} = 0.7 x (107 ~ 10%) A/m2 (2 8) 


If the capability or detecting the water flow exists, the 
direction of water flow can te known, and conseguently 
predict the epicenter of the earthquake. The earth current 
due to electrokinetic effects yields variations of the 
geomagnetic field on the surface. 

If the current flows uniformly in a horizontal 
circular cylinder of radius R, the magnetic field caused Dy 


the current is simply obtained [ky Biot and Savart's law, 
H = (I f 2724) °¢ 7{xrRe (2. 14) 


where r is the distance from the center of the cylinder to 


the point at which the magnetic field is observed. If 

ES eee (2. 15) 
then the field strength becomes 

H = 0.2 ~ 200 ¥ (20 


where the field strength is independent of R and L for this 
assumed geometry. Interestingly enough most observed magni- 
tudes of magnetic anomalies associated with earthguakes are 


almost all in the range estimated above. 


C. SOURCES OF ELECTROMAGNETIC FRECURSOR EVENTS 


1. Electric Fieids and Atmcspheric Electric Potential 


At the earthguake moment and just before it the 


effects due to disturbances of natural electromagnetic earth 


meera (NEMZE) were observed repeatedly. These effects 
involve quasi-stationary disturbances of the atmospheric 
electric potential (AEP) up to 1000 V/m, at distances up to 
some hundred kilometers from the epicenter. In the course of 
field observations the AEP perturbations were noted some 
hours before the earthguake moment within one or two hours 
period [Ref. 11]. 


AEP disturbances estimated for the charged surface 


are: 

Eaep =O 7 2 €vE,, (Zt) 
where 

Ge = ¢g N S (2.18) 


Under the assumptions that the linear charge density q, the 


surface density N, and the dislccation path length are: 
gq = 10°14 - 10719 C/n (22) 
N = 109 - 1022 ne (2n20)) 


AEP disturbances of 100 - 1000 V/m can be expected in the 


final stage of the earthquake precursor period. 
2. Magnetic Fields and Animal Behavoir 


The amplitude and time duration of the most charac- 
teristic time varying geomagnetic disturbances occurring at 
the earth's surface are listed in Table I . Another way of 
presenting the same data is in a plot of the spectrum of the 
total magnetic field. Such a spectrum is shown in Figure 
Zee Rete: 12), where skin defth (Z = 500 (.f£)7~% 5) is 


plotted as a function of frequency £ and conductivityc. It 
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TABLE 1 


Description of 
Phenomena 


Sub-audio frequency 
Fluctuaticns 


Micro-pulsations 
Giant pulsations 
Sudden ccmmencements 
Crochets 
Bay Phenomena 
Diurnal variation 
Storm Tine effects 


Twenty-seven day 
recurrence effects 


Annual variation 


Solar cycle secular 
variations 


bad se 


is evident that 


semi-diurnal, 


sun) and high frequency changes that are, 
tationary, that is net frequency stable. 


represent also a crude measure 


detection by 


animals can somehow understand 
well below 


Schemes for doing this are 


to detect changes 


ments presently in operation on 


twenty-seven day (due to 


ahimals. 


the most dominant signals 


it 25,5) 0m 


course, possible 


the background 


Time Variations in the Geomagnetic Field 


Typical 
Magnitude 
(gamma) 
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are the diurnal, 
the rotation of the 

in general, nons- 
This spectrum could 
of discrimination limits for 
that 


these variations well enough 


variations. 


necessary also for field experi- 


the San Andreas Fault where 


information on the large scale of these distrubances is used 


to enhance detection 


capability for 
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Figure 2.12 
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VAELAtLONS = Ret toy If animals can detect effects cecur— 
ring days before an earthquake, it is possible that these 
effects involve systematic change which would enhance the 
spectrum over a critical frequency range. 

There are a number of physical mechanisms that mght 


be considered as sources for earthquake-related animal 


| TABLE Il 


Order of Magnitude Estimates of the Various Possible 
Physical Mechanisms for Earthquake-Related Field 
Anomalies 
| A. MAGNETO HYDRODYNAMIC 
| ADBrnsecen = BOTwvV d > 0.1 gamma 
(o- = conductivit = 4 x 1077, v = velocit 
| D = length Scale} oo : ‘[ 
Be. TELLUBEC CURRENTS 
| ABrNAvcER = OE / D ~ 0.01 gamma 
(7 = conductivity, E = gecelectric field, 
| D = depth of current flow) 
C. PIEZO MAGNETIC EFFECT 
| AB= 40 S I~ 1 gamma 
| 


(S = stress change, I= magnetization) 


D. LOCALIZED CRACKS AND FRACTURES 
[oe = 


Oe gee ee ee came Os cers enrages cee UO pls a asec See eR pee ee ny Sem CE ey Ge, ee 





behavior. Table II lists several of these together with 
order of magnitude estimates of their amplitudes. It would 
appear that for these effects at least, it is difficult to 
generate source amplitudes of more than afew tens of 


gammas. The last process suggested which involves this 
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concerns the relatively unknown electrical and magnetic 
effects of local surface cracks and fractures might bea 


possible exception to this. 


Sa) 


Tift. ELECTROMAGNETIC PROPOGATICN CRITERIA FOR EARTHQUAKES 

A geological fault line becomes a source of acoustic, 
thermal and electromagnetic radiations when the tectonic 
plates move relative to each cther. The region near the 
fault line, shown in Figure 3.1 1S in both shear and 
compressive stress and the relative permittivity, Ep, , 
greater here than in the unstressed region having permit- 
tivity Er, - The fault width of 2Yo is the region of uniforn 
stress and the regions beyond this point are unstressed. 

When fault line motion commences, a voltage is generated 
causing current flow across an interface between two 
different materials. The interface can be formed by the 
following means: 

a. Forces near 100 Kbars can convert the insulating rock 
to a semiconductor {or conductor) in solid form [Ref. 14] so 
that charges can flow at the insulating-semiconductor inter- 
face which is also the voltage generating mechanism Since 
the work functions differ in the two media. 

b. Forces above 100 Kbars may cause plastic rock flcw to 
produce a tJliguid-solid interface [Ref. 15] with differing 
work functions causing current flow and charge storage in 
the solid portion having insulating properties. 

c. Shock waves can distort the dipole fields in the rock 
to produce a polarization in the material and increase the 
permittivity by nearly an order of magnitude [Ref. 16]. The 
transient charge is stored on the lower permittivity regions 


and accumulates with subsequent shocks. 
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Figure 3.1 Cross-Section of Earthquake Pault Line 
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Ae THE ELECTRIC FIEDDSESue eae 


The source of the electrcmagnetic radiations is the 
discharge of the accumulated electrons that were formed by 
shock waves or increasing pressures on the fault. A contin- 
uous corona discharge can occur from the nonuniform fields 
that form on Sharp points of the rock or a sudden discharge 
can result when the field exceeds the breakdown strength of 
the rocks in uniform or nonunifcrm fields. If the breakdown 
occurs in microseconds, radio waves in the megahertz region 
are radiated. This increasing pressure at a fault can be 
related to the strength and freguency spectrum of radiations 
which can be correlated to the acoustic data presently being 
recorded at the earthguake sites. 

The fault hypocenter, shown in Figure 3.1 is usually 
below the surface and is the scurce of the electromagnetic 
radiation, shown aS a point source. The emissions are frob- 
ably elliptically polarized but only the parallel field is 
depicted for clarity. Figure 3.2 shows the electric field 
vector, E in volts/meter in the incident X-Y plane and the 
magnetic field vector, H in apperes/meter, normal to the 
incident plane. The energy propcgates on the transverse mode 
(TEM) and the Poynting vector, Pin watts/square meter, on 
the X axis. 

The fault eficenter is shown in Figure 3.1 and is the 
location of the hypocenter at the fault surface where the 
electromagnetic radiations are a maximun. 

The electromagnetic radiaticns can be obtained from the 
retarded potential method but the objective is to only indi- 
cate the attenuation of the field with distance which can be 


found from the Maxwellian equations: 


ou =m 

VxexH=DE+E JE/Soat (3.1) 
aX 

Vx E=-poH / dt (3.2) 
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Figure 3.2 


VoLTS /METER 


Electromagnet Vectors at Fault Hypocenter 
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The parallel E vector mode, where the E field is in the 
incident plane, can be obtained from the vector equations to 


obtain the scalar differential equations: 
- QHZ / 9x = oEY + € SEY /at (3.3) 
aky /eoe——— fe Az 7 Xe (3.4) 


Let the arc discharge have the spectrum f1< £ < £2 and 
JWwt 


cs= 297 ff. If E = Eoe e , the scalar field reduces to: 
aEy / AX - y? Ey (x,t) = 0 (3.5) 
where 
yo Hyer” we “- 
Assume a scalar solution with E = Eoe eX ° Qe 2k | 


This gives the phase shift and attenuation: 


i (pw) 5 (jo /we -1) (3.7) 
=Flpw)'S (1- jo Aue)’ (3.8) 
= jlo) S(1-jo7 2wet «..) (3-9) 
= jlyw) S(1-j O/ 24) (3. 10) 
= jie) 8 + O/2 we (3. 1H 


The radiating field has the forn: 


Se 7 
~4 Jive 2 XjLPto . Ace 


E = Eo e (3.12) 
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Be. SKIN DEPTH AND RADIATED FRECUENCIES 


Mm Figure 3.1, the fault line region which has a 
constant compressive and shear stress is region one, and the 
unstressed region, away from fault line, designated region 
two. The permittivities and conductivities in the two 


regions are: 


oa aCe (3.13) 
wD, (3.14) 
é,= Ene. (3. 15) 
pes oe (3. 16) 


Where € is the permittivity of free space and Ca andé,, 
Meemene relative permittivities. The skin depth, $, is the 
distance where the radiation is reduced by the factor of Te 


from the initial value. The skin depth can be obtained froa 


equations: 
9 = 246, /07 = 206,40 (er) 
Using: 
f = 5x107 Hz, “= 100% x 10 Hz, BD; =10-S, (Geng) 
Cr = Se Gp = 8.8 aM0712 (3.19) 
S = 13.9 Kn (a2on 


Since the San Andreas fault is about 15 Km deep, the 


Skin depth is less than this value which indicates that 
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under the worst case of the hypccenter at the extreme depth, 
the electrical field is reduced slightly more than 1/e, but 
easily measured at the epicenter. 

The values of permittivity and conductivity vary with 
frequency and the upper limit c£ propogating frequency can 
be eStimated from the graph by Kraus and Carver [Ref. 17]. 
Using the worst case of urban ground at 100 megahertz, the 
skin depth 1s about 0.8 Km and if the hypocenter is halfway 
down the 15 Km tectonic plates, the attenuation is e¢6785, 
This limits the radiation to the range 5 X 107 < £ < 108 Hz. 


C. THE PROPOGATION CRITERION 


Figure 3.1 shows the radiation source at the hypocenter 
and the electric field, parallel to the incident flane, 
makes an incident angle 91 with the normal to the fault. The 
reflected wave 1S at an angle Or with the normal. The 
stressed region has a width 2Yo and a ae If the 
regions are insulating and lossless, the reflection coeffi- 


cient is given by: 


C= ar 7 Es (3-215 


wy oe} - | 2- sivta; 


- é ; (3. 2p 
COS &, 1h é, a Sin * 6) 

The reflected wave, Er, eguals the incident wave, Ei, 
when the radical in the numerator vanishes. The critica 


value of the incident angle for complete reflection is: 
singgics—"(227<)) => eee ae (3 235 


This relation obtains for either perpendicular or 
parallel polarization of the electric field and for all € 2 
@ic. Using Sia 25.0 for the stressed region and €r, = 5.0 


for the unstressed portion, gives 


qu 


81 = Sin-1(0.2)-5 = sin~1(0.45)- 26° (3. 24) 


Since the critical angle is small, the majority cr the 
radio emissions are reflected from the unstressed region for 
hypocenters deep in the fault. This wide region of radio 
emission may account for the abnormal human and aninal 
behavior distant from the fault line but within the stressed 
fault region. When the field in the air space above the 
epicenter exceeds 30 Kv/cn , the aixr breaks down and 
confirms the reports of lightning before and during large 


earthquakes. 
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IV. DESCRIPTION OF RESEARCH 


A. DESIGN GOALS 


One of the most important problems associated with using 
electromagnetic radiation aS a frecursor of earth quakes is 
how to select the point of observation. Placement of the 
system far from industrial areas, but near potential earth- 
guake activity was desired. A preliminary field survey of 
potential site locations was ccnducted under the direction 
of the United States Department of the Interior, Geological 
Survey. Site location along the San Andreas fault was chosen 
because of it being known historically and recently as an 
active earthquake area. The chcsen observation site is 1000 
feet from the San Andreas Rift Zone (fault) located at 121° 
23.5' W and 36° 45.5' N which is three miles east of San. 
Juan Bautista and seven miles scuth of Hollister. 

It became very important to secure spot window frecguen- 
cies for the system so not to have any interference fron 
transmitted radio waves. The Amy Frequency Coordinator for 
the Western United States identified and assigned to this 
project VHF bands “(38e85.4HZ Sande 1 50— 75 uma These 
frequencies, with a 50 KHz bandwidth, have been allocated 
With exclusive use for a period of one year. 

For observation of Signal intensity and direction in the 
VHF range, a data recording system was positioned with its 
directional sensors (antennas) looking both up and down the 
fault as shown in Figure 4.1 . The antennas (eight Yagi's) 
detect both horizontally and vertically polarized emissions 
at both reserved frequencies. | The omnidirectional antennas 
(two) are erected for detecticn at both freguencies, and 
with vertical polarization. Fach of the ten antennas are 


connected to a detector (communications receiver). 
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polarized and the other is horizontally polarized. 
There are two omni antennas vertically polarized at C. 
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Figure 4.1 Position of Antennas Along the San Andreas Fault 
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Figure 4.2 Block Diagram cf Data Reception Systen ~ 
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Low cost modularized receivers were procured and modi- 
fied to receive at the assigned frequencies. The output of 
the receiver is offset to nullify any ambient noise level 
and is fed into slow speed chart recorder to collect any 
detected emissions as a function of time. A block diagram of 
the measurement system is illustrated in Figure 4.2 . The 
noise measurement system can operate from 12 volt DC storage 
batteries or a power supply. The chart paper and storage 
batteries have to be replaced periodically. Electromagnetic 
emission data collected through this system will be related 
to seismic guake data provided by the United States 


Geological Survey Office, Palo Alto, Ca. 


Be SYSTEM COMPONENTS 


1. Receivers 


The VHF Converter Module 1s designed to amplify and 
convert the specific frequency to the HF range. The 
converter is supplied with a 10.7 MHz intermediate freguency 
(IF) transformer for the mixer output circuit. The mixer 
output circuit consists of a slug tuned coil, a capacitive 
voltage divider, and an RCA jack to provide a 50 ohm output. 
Figures 4.3,4.4 show schematic diagrams of the VHF Converter 
Modules. 

The IF Audio Module comprises a highly sensitive and 
selective IF amplifier, AM detector, audio amplifier, and 
sgueich system for use with the converter to make a VHF 
communications receiver. The selectivity of the unit was 
tailored to operate in the noise measurement system. The 
unit has automatic gain control (AGC) Omeooth LF and RF, 
with the RF AGC being delayed fcr best response. An integral 
log detector provides drive fcr the "S" meter (signal to 


Rustrak recorder) outrut. 
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VHF Converter Schematic for 38.45 MHz 


Figure 4.3 
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VHF Converter Schematic for 150.75 MHz 
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Figure 4.4 


The 10.7 MHZ Signal 1S converted to 4555582 “ang 
amplified by three stages. Each stage consists of several 
transistors biased from an integral zener regulated source. 
An optimally coupled three stage IF filter provides selec- 
tivity against adjacent channel interference. Figure 4.5 
Shows a schematic diagram of the IF Audio Module. Appendix 
A contains Table IV and Table V, data of receiver signal 
output voltages for a given RF signal input voltage at both 
VHF assigned frequencies. The output capabilities of the 
receiver were determined from this information on the 


S-meter (output signal to the chart recorder). 


2- Recorders 


The Rustrak DC recorder prints through the impinging 
action of its stylus driven ry the chopper bar against 
pressure-sensitive chart paper. Its presentation 1S a 
series of dots appearing aS a continuous line. Writing 
speed varies with motor speed. Chart speed and trace density 
depend on the ratio of the interchangeable gear box which 
couples the paper drive to the motor. Table III shows the 
chart speed and trace density. Figure 4.6 is a schematic o£ 
the Model 2146 dual trace chart recorder and Figure 4.7 isa 
schematic of the same nodel with a Signal processor included 
as shown. Both versions of Model 2146 had to be modified for 
Single channel use. The single channel modification (S-MOD) 
involved bypassing a switch (S-2) by putting a small jumper 
wire between the normally closed (N/C) and the normally open 
(N/O) terminals as shown in Figure 4.8 . Three of the 
recorders have signal processors to desensitize the signal 
input to one volt per division. The desensitizer board was 
demodified by rearranging the bias stepper connecticns at 
the Jones plug making non-bias stepper connections as shown 
in Figure 4.9 . There were no modifications necessary for 
the Model 388 chart recorder. 
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IF Audio Module Schematic 


Figure 4.5 
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TABLE IIIf 
Chart Speed and ‘Trace Density 


oz 
| 


Model Drive Writing Z 6 Gear Train No. 
Motor Speed S07 Oe Gear Ratio 
2194 1 1 strike/ 3 in/hr Chart Speed 
4 seconds 76 mm/hr Chart Speed. 
252 hrs Paper Duration 





(RPM) 
| 


388 2 1 strike/ 2 in/hr Chart Speed 
2 seconds 50 mayhr Chart Speed. 

378 hrs Paper Duration 
Trace Density 900 300 Strikes/iInch 
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3. Interface Amplifiers 


By design the receiver cutput voltages were approxi- 
Mately between 280 nV to 320 nV with RF inputs from approxi- 
mately 10 uV to 100 uv. Ideally this 40 mV range would drive 
the recorder from zero to full scale, but because of an op 
amp voltage offset only 80 % scale readings were achieved. 
An interface amplifier was designed to set full scale 
Minimum and maximum needle movement in the chart recorder. 
Calibration of the meter movement was accomplished by vari- 
able resistor settings as shown in the schematic of Figure 
4a l One 


4. Timer 


To compare seismological disturbance data with the 


noise measurement site recordings accurate time needed to be 


54 


| TU td 





H™ 
aD 
AID rn 
ON, 
a eee eee Biot v 7 
wna *T39 / AM 
2 
/ MM Sa1g saveL 
Ta 2 
Jaz ya 
gaa/im iv 
YIN Sy 





O) O 


DUA TI 
SSS WVIGVID 


aS, SS 
— er te EES ES AD a AEP a te ES AD aD a aS aa Se ST ST LT a ae 


Rustrak Mcdel 2146 Schematic 


Figure 4.6 
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Rustrak Model 2146 with Signal Processor Schematic 


Figure 4.7 
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Figure 4.8 Single Channel Modification (S-#OD) 
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Figure 4.9 Bias Stepper 
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Interface Amplifier with Voltage Generator 


Figure 4.10 
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indicated on the chart paper. Because of a requirement for 
ultra-low dissipation and operation using a battery source, 
a CMOS 21 stage counter timer (RCA CD4045A) was used in 
conjunction with CMOS 8 stage fresettable synchronous down 
counters (RCA CD40103B) anda CMOS dual J-K master-slave 
flip-flop (RCA CD4027B) to obtain the time. The 2.097152 MHz 
crystal establishes the one second clock pulse as shown in 
Figure 4.11. The time and the length of the pulse duration 
are set using 8-bit down counters. Waoen the timer pulses, 
the length of the pulse indicates the time. A fifteen minute 
pulse will last for six seconds. A one hour pulse will last 
for twelve seconds and a twenty-four hour pulse will last 
for eighteen seconds. Time can be initialized as necessary 


by the reset switch. 
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Time Tick Generator 


Figure 4.11 
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V- PERFORMANCE OF SYSTEM 

The informational parameter of a possible seismic event 
is variations in the level of the electromagnetic radiation 
envelope or in the number of radiation impulses whose ampli- 
tude exceeds a pre-chosen discrimination threshold. The 
noise measurement system is deSigned with a sensitivity 
range of 0.1 uV to 100 uV RF input. The system was given a 
"Shakedown test" on the roof of Spanagel Hall at the Naval 
Postgraduate School to see if it was sensitive to variations 
in the level of electromagnetic radiation based on a 24 hour 
period. 

Along the San Andreas Fault, it is expected that varia- 
tions in Signal level will occur due to seasonal chanyes. 
For instance, during autumn and winter the highest signal 
level falls into the sunless part of the day. Springe 
summer, andthe beginning of autumn have two maxima during 
the 24 hour period, the nightly and the post-noon ones. 
Figures 5.1, 5.2, 5.3, 5.4, and 5.5 are recordings wii 


Show the noise levels over a 24 hour period. 
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Signal Level Data from Noise Measurement Systen 
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Fagumes 5. 1 
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Figure 5.2 
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Figure 5.3 
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Figure 5.4 
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Figure 5.5 


67 


VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


The design, construction, and testing of an RF noise 
measurement and recording system in the 30 MHz and 150 MHz 
range has been described. Test results show that background 
noise levels can be received and recorded from tenths of 
microvolts to volts. Theory has shown that anomalous elec- 
tromagnetic radiation levels from earthquakes can occur. 
The noise measurement system data will not only detect the 
presence of increased background noise levels at the desig- 
nated frequencies, but it may also give an indication of the 
direction of reception relative to the San Andreas Fault due 
to the polarization of the antennas. | 

Besides physically erecting the antennas along the 
fault, the most difficult problem was designing a receiver/ 
recorder system that would be sensitive to slight changes in 
background noise level at frequencies near those of citizen 
band and amateur radio operators. The signal output of the 
receiver was sent through an interface amplifier and the 
Chart recorders modified to acccmplish this. 

Appendix B presents electromagnetic emissions and other 
observations possibly related to earthquake precusor 
activity. Measurements of electromagnetic radiation from 
earthquakes were detected by changes in background noise 
levels, but in the very low, low, and high frequency ranges. 
The results of this research suggest that electromagnetic 
radiation measurements has pctential as an earthguake 


prediction tool. 
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B. RECOMAENDATIONS 


Data acquired from the VEF noise measurement systen 
should be taken for at least a year observing background 
noise level variations. Since seismological data can be 
acquired from the Geological Survey Office, any Significant 
earthguake in the area of the noise meaSurement Site should 
be immediately compared with the recorded data. Otherwise 
data can be examined on a periodic basis. 

Although several studieS were presented concerning the 
mechanism by which rocks emit radiowave energy when they are 
crushed, further investigation is needed in this area and 
the following: 

1. The relation between the intensity of the background 
hoise inthe VHF range and the location of the 
epicenter. The antennas of the noise measurement 
system are looking up and down the fault which makes 
location protable. 

2. The relation between the intensity of the precursor 
VHF radio noise and the earthquake's magnitude and 
depth of focus. The data obtained may show similar 
variations in background noise levels at 38 MHz or 
150 MHz or both for earthguakes of similar magnitude 
and depth of focus. 

3. Long term data collection is needed to  oftain 
through order-of-magnitude calculations, the amount 


of VHF radio noise radiated in these events. 
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RECEIVER SENSITIVITY DATA AT VHF FREQUENCIES 
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RF Input in av an Signal Output in av 








RF S RF S RF S RF- S 
fn fos Toss Yoon foe ese 
Cra (0.226 2 26 Om ) 0.298 O. p0.324 350 fo. 0.344 | 
OF 6 


=e [es [oa [ef 
seats pepe ee oe 


Otte, We ee, en 


_ 


APPENDIX B 
ELECTROMAGNETIC EMISSICNS FROM EARTHQUAKES 


A. COALINGA 


An ionospheric perturbation that was produced by the 
Coalinga earthguake of 2 May, 1983, was detected by a 
network On. High frequency radio links in northern 
California. The ionospheric refraction regions of all five 
HF propagation paths, at distances between 160 and 285 kn 
(horizontal range) from the epicenter, as Shown in Figure 
B.1, were all affected by a ground-motion-induced acoustic 
pulse, aS Shown in Figure B.2, that propagated to ionos- 
pheric heights. These observations appear to be the first 
ionospheric disturbances to be reported this close to an 


earthguake epicenter [Ref. 18]. 


B. RUSSIA 


Figure B.3 illustrates one example of the anomalous 
increase in electromagnetic radiation possibly related to 
earthquake activity on September 16, 1978. In this case, the 
epicenter of earthguake was located in central Iran, and the 
receiving equipment was set up in a tunnel 50 m_ below the 
ground surface at Caucasus, USSR, 1200 Km distant. The 
on-set time of electromagnetic radiation at frequencies of 
27 KHz and 1.63 MHZ was about 30 minutes before the main 
shock with magnitude of 7.4, as shown in Figure B.3 
Hert. 2}. 
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Figure B.1 
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Figure B.3 Electromagnetic Noise Levels at Caucasus 
C. JAPAN 


On March 31, 1980, an earthquake occurred at 1633 Japan 
Standard Time (JST). The magnitude was 7, and the epicenter 
was located at 35.49 N and 135.39 E in the Kyoto prefecture. 
The focus of this earthquake was very deep, at approximately 
u80 Km. The distance between the location of the station and 
epicenter was around 250 Kn. 

As shown in Figure B.4 {Ref. 2], about half an hour before 
the shock, the 81 KHz data reception system recorded an 
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Figure B.4 Electromagnetic Noise Levels at Sugadaira 


anomalous increase in the intensity. The amplitude exceeded 
the normal level by more than 15 dB. The normal level was 
approximately 8 u¥V. It can also be seen in Figure B.4 that 
the intensity dropped sharply back to the previous level 
exactly at the moment of the shecck. 

Another event occurred at 1247 JST on January 28, 1981, as 
Shown in Figure B.5 [Ref. 2]. ‘The location of the epicenter 
was in the western Ibaragi prefecture on the northeastern 


side of Tokyo, and the distance between the epicenter and 
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Figure B.5 Electromagnetic Noise Levels at Suginani 


the observation point was approximately 50 Km. The depth of 
focus was 60 Km and the magnitude was 5.0. The noise lovel 
at 81 KHz increased to about 12 dB above the threshold level 
about 45 minutes before the main shock. Such a feature is 


quite Similar to the typical events presented in Figure B.4,. 
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